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Abstract 

Three elementary canonical transformations are shown both to have quantum 
implementations as finite transformations and to generate, classically and infinites- 
imally, the full canonical algebra. A general canonical transformation can, in prin- 
ciple, be realized quantum mechanically as a product of these transformations. It is 
found that the intertwining of two super-Hamiltonians is equivalent to there being 
a canonical transformation between them. A consequence is that the procedure for 
solving a differential equation can be viewed as a sequence of elementary canonical 
transformations trivializing the super-Hamiltonian associated to the equation. It is 
proposed that the quantum integrability of a system is equivalent to the existence 
of such a sequence. 
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Canonical transformations are a powerful tool of classical mechanics whose strength 
has not been fully realized in quantum mechanics. Despite work by Dirac[|l| and others]^, 
^, ^, ^, there is a persistent distrust of the utility of canonical transformations in quan- 
tum mechanics. The general feeling lies between believing they can't be used in quantum 
mechanics and thinking that at best they are of only limited interest. The truth, how- 
ever, is that effectively all exactly solvable problems in quantum mechanics are, or can 
be, solved by canonical transformations. Each of the familiar tools for integrating a dif- 
ferential equation — change of variables, extracting a factor from the dependent variable, 
Fourier transform, etc. — has a realization as a canonical transformation. Only procedures 
whose origins lie in approximation, such as power series expansion, are not essentially 
canonical transformations. 

By making the quantum implementation of canonical transformations explicit, the 
mechanics of solving differential equations is made more transparent. This is especially 
useful in problems of more than one variable where notational complications begin to 
cloud insight. Furthermore, the usefulness of canonical transformations is not limited 
to exactly solvable problems; they can also be important in preparing a problem for 
approximation. As well, by bringing to quantum mechanics a tool so central to our 
understanding of classical mechanics, our grasp of their relationship is deepened, and 
new lines of investigation are opened. 

Perhaps most significantly, canonical transformations provide a unifying structure 
in which to study quantum integrability. The method of intertwining 0, in which one 
constructs an operator D which transforms between two operators 

= DH^''^ (1) 

so that their eigenfunctions are related 

^a) = D^(o), (2) 

has already been proposed for this role0. This is because it underlies the factorization 
method of Infeld and Hullp, supersymmetry in quantum mechanics^, and Lie algebraic 
methods] T0[|, which together solve most integrable problems in quantum mechanics. The 



method of intertwining was limited in the past by the need to make an ansatz for the in- 
tertwining operator D. We shall see that the intertwining operator is simply the canonical 
transformation between 7i^^^ and Ti*-^-*. 

This observation considerably deepens one's intuition for constructing intertwin- 
ing operators. In addition, it explains why intertwining underlies the Lie algebraic 
methods] TI|]: the symmetry group preserving the form of a differential equation is a 



subgroup of the full canonical group. More, the infinite-dimensional Kac-Moody algebra 
associated, for example, with the soliton solutions of the KdV equation, is explained 
through the Lax formalism because it is a subalgebra of the full canonical algebra. 



2 



In the following, three elementary canonical transformations and their quantum im- 
plementations will be introduced. Using these, several additional composite elementary 
transformations will be constructed. It will be shown that, classically, the infinitesimal 
forms of these transformations generate the full algebra of canonical transformations. 
This implies that in principle we have the tools to make any quantum canonical trans- 
formation. 

In practice, we are limited to composing a certain collection of finite transforma- 
tions. These transformations will be recognized as familiar tools for solving differential 
equations: changes of variable, similarity (gauge) transformations, Fourier transform. In 
addition, the standard first order differential intertwining operator will be realized as a 
canonical transformation. It is noteworthy because it is a transformation whose action 
is simpler in the quantum context than in the classical. Finally, a few examples are done 
to illustrate various aspects of the use of canonical transformations. Many more compli- 
cated examples follow straightforwardly by use of the differential intertwining canonical 
transformation. 

As this paper was being prepared, the author became aware of recent work by Leyvraz 
and Seligman |]TB[ and Deenen which overlaps in part with the results given here. 



1 Formal aspects 
1.1 Introduction 

The canonical transformations discussed in this paper may be applied to any wave equa- 
tion. For convenience, attention will be focused on the non-relativistic Schrodinger equa- 
tion. The presumption is that the wavefunctions of an initial equation are of interest, 
and the intent is to construct them by canonically transforming to simpler equations. 

Since canonical transformations involving the time and its conjugate momentum 
Pq are useful, it is profitable to work in the extended phase space in which these are 



adjoined to the spatial coordinates and momenta|jT^. The Schrodinger operator is the 
super-Hamiltonian 

^(q,p) =Po + ^(%'9i'Pi)' (3) 

where q, p are used to denote all of the extended phase space coordinates. For notational 
convenience, only one spatial degree of freedom is indicated explicitly. The Schrodinger 
equation in the coordinate representation is the constraint equation 

^^(q) = 0. (4) 

A canonical transformation is performed by applying an operator, say P, on the left 
of the super-Hamiltonian constraint equation and rewriting the result as a new constraint 
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equation 

= P7^(0VW=7^(iV«. (5) 
This gives the transformed super-Hamiltonian 

7^«(q,p)=P7^(°)(q,p)p-\ (6) 

and the transformed wavefunction 

^(1) = p^{o)^ 

The operators P, P^^ can be brought inside the super-Hamiltonian to act directly 
on the phase space coordinates 

7i:«(q, p) = n^^\PqP-\ PpP-1). (8) 

This is possible for any function on phase space possessing a Laurent expansion about 
some point, as can be seen by expanding the function, applying the transformation and 
resumming. In practice it holds for all super-Hamiltonians of interest. Since the canonical 
commutation relations are preserved by this transformation 

[q,p] = [PqP-\PpP-']=t, (9) 

the effect is to transform the original phase space coordinates by the canonical transfor- 
mation, which was the object. Note that (^) is the same as the intertwining condition (|I|). 
This establishes the equivalence of the canonical transformation P and the intertwining 
operator D. 

The goal of constructing the wavefunctions of Ti.^^^ is attained if, by a sequence of 
canonical transformations, one can transform to an Ti.^^^ which is known to be integrable. 
It is conjectured that a super-Hamiltonian is integrable if and only if it can be canonically 
transformed to the trivial super-Hamiltonian 

n = Po- (10) 

The canonical transformation to a trivial super-Hamiltonian has special significance in 
classical mechanics: this transformation makes all the variables constants of the motion, 
and these constants may be chosen to be the initial conditions of the system. Inverting 
this transformation classically solves the equations of motion. 
The inner product is assumed to be conserved 

J dq = Jdq (11) 

(where q runs over a spatial hypersurface ||15|| ). This determines a transformed measure 
density, which in turn may be used to define the adjoint of P. Under this definition. 
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P^P = 1, and P is seen to be a unitary transformation. In the general case, the measure 
density may be operator-valued — consider, for example, the Wronskian operator appear- 
ing in the Klein-Gordon inner product. Details of the transformation of the measure 
density will not be discussed in this paper; they are straightforward to work out from 
the inner product. 

1.2 Infinitesimal Canonical Transformations 

Classically, an infinitesimal canonical transformation is generated by an (infinitesimal) 
generating function F(g,p)|[T^. Associated to this generating function is the Hamiltonian 
vector field 

VF = F^pdg - F^gdp (12) 
whose action on a function u on phase space is the infinitesimal transformation 

6fu = evpu = —e{F, u} (13) 

where {F,u} is the classical Poisson bracket. Through the correspondence between clas- 
sical and quantum theory in which Poisson brackets go over into commutator brackets 
times i, this canonical transformation can be expressed in terms of the quantum operator 
exp(— ieF) 

Because of operator ordering ambiguities in defining the quantum versions of F and u, 
the classical and quantum expressions for the infinitesimal transformation can differ by 
higher order terms in h (at the same order of e). 

In the equations that follow, quantum variables are not distinguished notationally 
from classical ones. Unless otherwise stated, all variables are handled according to their 
familiar quantum contexts, e.g. as operators in Hamiltonians and as c-numbers in wave- 
functions. Products of operators are assumed to be ordered as written; classical versions 
are obtained by allowing p and q to commute. 

There are three elementary canonical transformations in one- variable which have well- 
known implementations as finite quantum transformations. They are similarity (gauge) 
transformations, point canonical (coordinate) transformations and the interchange of 
coordinates and momenta. By finding the algebra generated by the infinitesimal versions 
of these transformations, the class of transformations that can be reached with them can 
be determined. The somewhat surprising result is that classically they generate the full 
canonical algebra. 

The significance of this classical result is that since each of the elementary and compos- 
ite elementary transformations has a quantum implementation as a finite transformation. 
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in principle any general canonical transformation can be implemented quantum mechan- 
ically. In practice at present, one is limited to finite products of the transformations. 
Nevertheless this is a very powerful tool for solving problems in quantum mechanics. 

Similarity (gauge) transformations are infinitesimally generated by Fs = f{q)- Point 
canonical transformations are infinitesimally generated by Fp — f{q)p. The discrete 
transformation I interchanging the role of coordinate and momentum is 

p = Ip'r^ = -q\ q = Iq'r^ = p. (15) 

Using the interchange operator, composite elementary transformations which are nonlin- 
ear in the momentum can be formed. They are the composite similarity transformation, 
infinitesimally generated by 

Fcs^IFsI-'^fip), 
and the composite point canonical transformation, infinitesimally generated by 

Fcp = IFpI-^ = -f{p)q. 

Each of these corresponds to a finite transformation through application of the inter- 
change operator to the finite forms of the similarity and point canonical transformations. 

The many- variable generalization of the similarity transformation is straightforward. 
It is infinitesimally generated by Fsn = /(?i, ■ ■ ■ , ?„)■ Interchanging any set of coordi- 
nates with their conjugate momenta gives the composite many- variable similarity trans- 
formations. For example, in two variables, one has the infinitesimal generating functions 

Fcs2a = f{Pi,q2) and Fcs2b = /(Pi,P2)- 

The observation is now made that classically the algebra generated by the elementary 
and composite elementary transformations is the full canonical algebra 

[vf,Vg] = ~V{KG}- (16) 

where F, G are arbitrary functions on phase space. It is to be expected that Hamilto- 
nian vector fields will have this algebraic structure. More surprising is that the above 
collection of generating functions produce a general function on phase space through the 
Poisson bracket operation. This is verified by taking commutators of the different types 
of transformations. 

Consider the two- variable case — the many- variable case follows similarly. Introducing 
the monomial generating functions 

F^L = ql^\2'^'pr'pr\ (17) 
where j, k,n,m e Z, the Poisson bracket of two of these is 

{^;^.,^^^J = ((i2 + l)(ni + l)-(ji + l)(n2 + l))^;S^t^?iVi (18) 



nimi ' n2m2 . 

+ + l)(mi + 1) - (fci + l){m.2 + l))F;Ji:a\ ttiL.- 
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Inspection shows that a general monomial can be constructed by beginning with the 
monomial forms Ti , F^Li F° z\, F~\^~^ , etc., which generate the elementary and 
composite elementary canonical transformations. 

By taking linear combinations, one can form any function having a Laurent expansion 
about some point (not necessarily the origin) . To avoid having to use formal Laurent ex- 
pansions about nonzero points for functions like g^/^ or In greater generality is obtained 
by working with generating functions of the form 

i^ = /l(?l)/2(?2)^(Pl)^(P2)- (19) 

This produces any function which can be represented as a sum of separable products. 



1.3 Quantum Implementations 

Each of the elementary canonical transformations can be implemented quantum me- 
chanically as a finite transformation. Their action is collected in Fig. 1, and each will be 
reviewed below. 

The interchange of coordinates and momenta 

p = Ip'r^ = -q', q = Iq'r^ = p'. (20) 



I - tAtt^. r die"^'" (21) 



is implemented through the Fourier transform operator 

1 

(2^ 
for which it is evident that 

Iq = -idqJ, Ip = -q'l. (22) 
The wavefunction is transformed 

^l;^'\q') = I^l;^'\q) = J^^ dqe"^'^^'\q). (23) 

The inverse interchange is 

p = q\ q = -p'- (24) 
It is implemented by the inverse Fourier transform 

I fOO 



The similarity transformation in one- variable is implemented by the operator 

Sm = e-^(^) (26) 
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where f{q) is an arbitrary function of the coordinates. This induces the canonical trans- 
formation 

P = Sf^q,)p'S~fl,^^p' -if^q,, (27) 
Q ^ Sf(^q,)q'Sf(^g,) = q'. 

The composite similarity transformation is implemented by applying the interchange 
operator to -S'j(q) to exchange coordinates for momentum in /. In the one-variable case, 
the composite similarity transformation operator is 

Sf(p) = ISf^,^I-'^e-f^\ (28) 

It produces the canonical transformation 

p = Sf(^,)p'S-fl,,^^p', (29) 

q = Sf^,)q'Sfl,^^q' + if^p,. 

In the many-variable case, the function / may involve either the coordinate or its 
conjugate momentum for each variable. Because the different variables commute amongst 
themselves, each responds to the (composite) similarity operator as if it were a one- 
variable operator in that variable, and the other variables are treated as parameters. 
Thus, for each coordinate (momentum) of which / is a function, the corresponding 
conjugate momentum (coordinate) is shifted as in the one- variable case. 

The finite point canonical transformation takes a different form than the infinitesimal 
version. The effect of the point canonical transformation Pf(q) is to implement the change 
of variables 

q = P;(,,)g'P;(;,) =/(?'), (30) 



/(</') - ■ 

J ,q' 

The effect of Pf[q) on the wavef unction is 

z^W(g)=P^(,)^(°)(g)=^(°)(/(g)). (31) 

The composite point canonical transformation is formed by composition with the 
interchange operator 

Pfip) = ^Pm^~^- (32) 

It has the effect of making a change of variables on the momentum 

q = ■^(^, (33) 

J,p' 

P = f{p')- 
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The operator ordering of the transformed q is determined by the action of the interchange 
operator on the coordinate point canonical transformation. 

A two- variable generalization is also of interest. It is defined by the canonical trans- 
formation 



qi = f{q'i,q2), (34) 
1 

z 



Pi = ^Pl 



,91 

^2 = q'2^ 

P2 = P2 + T^Pl- 

Using this, the result of exponentiating the infinitesimal generating function F = 
g{q2)h{q^)p-^ can be deduced. Since different variables commute, the action of e~*^ on 
is that of a similarity transformation (^) and we have 



Y^ = 9{q'2),,ih{q[)- (35) 

J .c 



Equating this with the transformation in (|3^) gives 

/, 

'<?i 

Defining 

the chain rule simplifies this equation to 

f,H - = 

which has the general solution 

f = f{H{q,)+g{q,)). 
Since the transformation must reduce to the identity when (7 = 0, we find 

f = H-\H{q,)+g{q2)). (36) 



This agrees with a result found independently by Deenen. 

As an explicit example, take g = const and h = q^. The infinitesimal transformation 
is then that of the Virasoro generators. The finite transformation is (for m 7^ 1) 

f = [ql-^- + [I - m)gf''-^\ (37) 
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while for m = 1 

/ = e^qi- (38) 

I remark that for m = 0, 1, 2 these transformations generate the group SL{2, C). 

In the one- variable case with g taken to be a constant, Eq. (|36D gives the finite 
transformation produced by the infinitesimal generating function F = gh{q)p. It can 
also be read as a functional equation 

H{f{q)) = H{q)+g (39) 

which is to be solved for H given /. Not surprisingly, this equation does not have a 
solution for all /. The implication is that not all point canonical transformations can be 
expressed as the exponential of an infinitesimal transformation. This is an explicit demon- 
stration of the well-known property of the diffeomorphism group that the exponential 
map does not cover a neighborhood of the identity[|r^]. This is a property of infinite- 
dimensional Lie groups and stands in contrast to the situation in finite-dimensional Lie 
groups. 

A corollary is that the product of the exponentials of two generators cannot always 
be expressed as an exponential of a third generator. For this reason, it is generally not 
useful to express point canonical transformations in exponential form, but rather to note 
directly the change of coordinate they produce. (As reassurance, it is true that every 
point canonical transformation can be expressed as a finite product of exponentials [0.) 



1.4 Linear Canonical Transformations 

The linear canonical transformations form a finite-dimensional subgroup of all canonical 
transformations, and there has been much interest in them in the context of coherent 
states [@, H, ^. As canonical transformations, they can be constructed from a product of 
elementary transformations. 

Consider the case of a single variable. A linear composite similarity transformation 

p = p", q = q"- - iap"-. 

transforms the wavefunction 

(A superscript is used to indicate the generation of the transformation. Roman letters are 
used to avoid confusion with powers of the variable. Subscripts are used to distinguish 
variables when necessary.) A linear similarity transformation 
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makes the change 

Finally a scaling of the coordinate 

7 

gives 

The full transformation is 

p= -p' - ip^q'^ q = p' + ^(1 - aP)q', (40) 

7 7 

with 

= e*i'^^^Pe-^^'/2e"^''/2^(°)(?)- (41) 

A general SL{2, C) = Sp{2, C) transformation is of the form p = ap' + bq', q = cp' + dq' 
where ad — be = 1. This gives the correspondence s = ic/a, [3 = iab, 7 = 1/a. 

By expressing ^p^^^ as the Fourier transform of 'ip^^\ an integral representation is found 
for ip^^^ which does not explicitly involve exponentials of differential operators 

= / rf^'em'.-/37^<?V2WV2^(o)^^/)_ (42) 

(27r)^/2 J 

A related result is given by Moshinsky^]. 

The operators p^, g^, {qp+pq)/2 generate a realization of the SL{2, C) algebra. Since 
SL{2, C) is a finite-dimensional Lie group, every element of the group can be expressed 
as an exponential of an element of the algebra. As well, a product of exponentials of 
elements of the algebra can be expressed as the exponential of another element of the 
algebra. Thus, a given linear canonical transformation may be expressed in many ways 
as a product of elementary transformations. Each will give an expression analogous to 

m or m- 

The generalization to many- variables is straightforward. The group of linear canonical 
transformations is Sp{2n,C), and a realization of it is found from the linear similarity, 
composite similarity and scaling transformations. Realizations of other finite-dimensional 
Lie groups in terms of canonical transformations are found by treating them as subgroups 
oiSp{2n,C). 

By expressing the coordinates and momenta in terms of harmonic oscillator creation 
and annihilation operators, one finds the expressions for the action of linear canonical 
transformations on coherent states^, These are useful for handling squeezed states 
in quantum optics^. 
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2 Applications 



2.1 Intertwining 

The most widely known example of intertwining is that of the Darboux transformation 
between two Hamiltonian operators in potential-form[^ 

H^'^ =p^ + V^iq). (43) 

Here, one would like to know what potentials can be reached from a given one by an 
intertwining transformation ([^). In the usual approach, an ansatz is made that the 
intertwining operator D is a differential operator, and the operator is constructed by 
requiring that the intertwining relation (|I|) be satisfied. It is satisfactory to begin with 
first order operators because the intertwining transformation can be iterated. 

Since the intertwining operator is a canonical transformation, it can be constructed 
as a sequence of elementary canonical transformations. The construction is interesting 
because it involves a transformation in which operator ordering makes the the quantum 
case simpler than the classical. This is the source of the power of the differential operator 
ansatz, and it is the first example of a transformation which favors the quantum problem 
over the classical. 

Begin with a Hamiltonian with potential 

H^^^ =p' + Vo{q). (44) 
The potential may be cancelled by making a one-variable shift 

p = p'' -igiq"), q = q''. 

This gives the Hamiltonian 
together with the Ricatti equation 

g,,a+g^ = V^ + \. (46) 

The transformed wavefunction is 

=e-/^"^«>(°). (47) 
A composite one-variable shift is made on the coordinate 

g» = g^-4, p^=p\ (48) 
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This is the key step, and it has such wide apphcation as to merit its own name, the 
differential intertwining transformation. The transformed coordinate can be expressed 
as 

This has the very interesting property 

9i^' - ^) = 9{pV^,) = 9{q') - (49) 

for functions with a Laurent expansion. Note that only the first term in the Taylor ex- 
pansion of g appears — classically the full Taylor expansion would have arisen. This is the 
essentially new feature of the differential intertwining transformation. The Hamiltonian 
is then 

i/W = - 2igp'' - 2g^gt - A. (50) 
and the transformed wavefunction is 

^W(g^) = e/*'/''V(")(g'') = -z9,.^(")(g''). (51) 

The one-variable shift 

= p' + tgiq'), q' = q\ 
cancels the term linear in the momentum giving the Hamiltonian 

i7(i)=/ + T4(g'), (52) 

with the new potential 

Vi = -g,,>+g'-X. (53) 

The transformed wavefunction is 

^W{q') = el'^'''^'^^'\q'). (54) 
In terms of the original wavefunction, this is 

^(i)(g) = eS9^^pe-S9dq^'^^)^q) = _ g)^(o\q). (55) 



Comparing (|^) and (|53D , the change in potential is 

V,-Vo = -2g^g. (56) 



The Ricatti equation (|^) [or (|53D] can be solved to find that g is given by the logarith- 
mic derivative of an eigenfunction of H^^^ (or the negative logarithmic derivative of an 
eigenfunction of H^^^) with eigenvalue A. This is the standard result from intertwining 0. 
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If one inverts ( ^5|) to obtain ip^^'' in terms of ip^'^\ an integral operator arises. To 
obtain a differential operator relation, one may note that taking g —g interchanges Vq 
and Vi- This implies that 

V;W(g) = -^(9,+^?)^«(g). (57) 

Alternatively, a different sequence of canonical transformations can be used which 
give an integral operator relating ip^'^^ to ip^^'^ which becomes a differential operator upon 
inversion. Beginning from (ID, the one- variable shift p = p"" + igiq""), q = q"' gives 
the Hamiltonian 

where Vq satisfies (|^). Note that has been ordered on the left of g. This is to facilitate 
the transformation 

g" = g' + 4, P' = p'. 
p° 

The transformed coordinate has the property that 

i _ 1 

This leads to the transformed Hamiltonian 

ff{b) ^ ^ 2ip''g - 2g^gt - A. (60) 



q' + - = -qV. (59) 



A similarity transformation p^ = p' — ig{q'), q^ = q' cancels the linear momentum term, 
leaving the final Hamiltonian 

i7«=/ + \4, (61) 
where Vi is given by (]5^). The final wavefunction in terms of the original is 

^(1) = e-/^'^V'e/^'^^^(°n?)- (62) 

Inverting this gives the expected differential relation (^). 

The method of intertwining has been realized in terms of canonical transformations. 
This means that all of the problems which can be solved by intertwining can be solved 
with canonical transformations. In particular, this means all problems which are essen- 
tially hypergeometric, confluent hypergeometric, or one of their generalizations. 
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2.2 Non-relativistic Free Particle 

The non-relativistic free particle is an easily solved problem which doesn't require any 
sophisticated machinery. It may however serve to illustrate a number of features of 
the use of canonical transformations, and, for this reason, it will be treated somewhat 
exhaustively. 

The free particle super-Hamiltonian is 

7i(o)=PQ+p2 (63) 

It may be immediately trivialized 

7^^"^ = pI (64) 
by the two- variable similarity transformation 

Po = - P"' P = P" 

?o = ?o ? = ?" + 2p"?^ 

The original wavefunction is given in terms of the transformed one by 

^(o)(g,g^j) = e-^f'«o^('^Hg). (65) 

where the solution of the trivialized super-Hamiltonian constraint TC^^^^l>^°'^ = is any 
^Q-independent function. This formula is just the formal expression for the evolution of 
an initial wavefunction in terms of the exponential of the Hamiltonian. In general, this 
formal result is insufficiently explicit, but for the free particle, it can be used to find more 
useful forms of the wavefunction. 

For example, if the initial wavefunction is taken to be a plane wave ■0^"^ = exp{ikq), 
one finds the plane wave stationary solution 

^(0) ^ ^ikq-ik\_ (gg^ 

If the initial wavefunction is a delta function at x, then using the Fourier integral repre- 
sentation of the delta function, the wavefunction is 

27r J—oo 

Acting with the operator inside the integral and integrating the resulting Gaussian gives 

= (47rzgo)"'/'e'(''-^)'/^''o. (67) 
This is the free particle Green's function. 
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It is clear that the nature of the wavefunction depends on the initial wavefunction 
used to generate it. This obvious comment is important to bear in mind because the 
"natural" solution of a transformed problem will not always correspond to the desired 
solution. To see this, consider a second approach to the free particle. An interchange 
t r ansf ormation 

p = -q" 

q = p", 

is equivalent to taking the Fourier transform of the original super-Hamiltonian and gives 

H^-) = Po + g"'- (68) 
The original wavefunction is given by the inverse Fourier transform 

^^°^(^) = j_Jqe-'''^^^'^{q''). (69) 

The new equation H^'^^^p^''^ = can be solved "naturally" in two ways. The first is 
simply to integrate with respect to ^q. This gives 

where /(g") is an arbitrary gg-independent function that arises as an integration constant. 
One finds the wavefunction 

1 poo 2 

= -(^2 dq^e-^'^'^''-'^'' "ofiq-). (71) 

This is a less familiar form of the wavefunction — it is of course just a momentum space 
representation — and it is perhaps not immediately obvious how to obtain the solutions 
above. Inspection shows that if / = 5(g" + k), one finds the plane wave stationary 
solution. Alternatively, if / is taken to be 

then, at = 0, one has ip^'^^ = S{q — x) and, evaluating the integral, one finds again the 
Green's function (|67D. 

The second "natural" approach is to separate variables 

=0(«)(g'^)e-^'='''o. (73) 
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This results in the equation 

(g"' - = 0, (74) 

which has as its solution 

0W = - k) (75) 

[k can have either sign). Now inverting the interchange operation gives the familiar plane 
wave solution 

= I ikg-ik^qg (jQ\ 

y (27r)i/2'' ■ ^""^ 

The Green's function solution is however no longer immediate. 

The first of these approaches can itself be implemented as a canonical transformation. 
The similarity transformation 

P^o=Po~ f p'^ = / - 

go = 4 <f = g' 

trivializes the super-Hamiltonian 

n^'^ = pI 

and gives the wavefunction 

^(a) _ e-g'^''?S^W. (77) 

The wavefunction ^^''^ is any gQ-independent function. The result for the original wave- 
function is then ([711). 

The second approach of separation of variables is not so much a canonical transfor- 
mation as a realization of the assertion that the solution space of the super-Hamiltonian 
has a product structure. 



2.3 Harmonic Oscillator 

The harmonic oscillator is the paradigmatic problem in quantum mechanics, and it is a 
test piece for any method. Its solution by canonical transformation reemphasizes how the 
form of a solution is affected by the details of evaluating the product of operators repre- 
senting the canonical transformation. Also, it is observed that more than one canonical 
transformation to triviality is needed to obtain both independent solutions of the original 
Hamiltonian. 

The Hamiltonian for the harmonic oscillator is 

ff(o)=p2^^Y- (78) 
A one-variable similarity transformation 

p = p" + ituq", q = q" 
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will cancel the quadratic term in the coordinate leaving 

This is recognized as the equation for the Hermite polynomials. However, since solution 
by power series expansion is not a canonical transformation, so much as a method of 
approximation, additional transformations are needed. The composite one-variable shift 

cancels the quadratic term in the momentum, giving the Hamiltonian 

H^^^ ^2iujqy + (79) 
Finally, the point canonical transformation 

= e-'^'p', q^ = e«'. 
transforms this to action- angle form 

H^^^ ^2iu;p' + 00. (80) 
In terms of the super-Hamiltonian, this is 

7^(1) = p(, + 2iu;p' + u. (81) 
A final two- variable similarity transformation 

P'o = Po - 2«^P" - (l'o = Qo 

p'^p", q' ^q'^ + 2iu;q'^ 

trivializes the super-Hamiltonian 

n^''^ = pI (82) 

The wavefunction ip^^^ is given in terms of by 

i^^'\q) = Pe.S,2/,^S_^^2/,i;^'\q), 
which may be inverted to find 

^^^H?) = S^,y2S-j,y,u^Pin,^^'l (83) 
Prom the eigenfunctions of H^^^ 

^(1) _ ^nq-i{2n+l)uiqQ ^ 
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one has 

This is the correct (unnormaUzed) harmonic osciUator eigenfunction. This formula is 
vahd for complex n. Requiring that the wavefunction be normalizable fixes n to be a 
non- negative real integer. For other n, one finds an infinite power series in 1/g which is 
divergent at g = 0. 

As remarked above, i?^") is the Hamiltonian whose solutions are the Hermite poly- 
nomials. This implies 

HM oc e-^^^^'V (85) 
Given this form, it is immediate that is the lowering operator 

dqHniq) (xnHn_i{q). (86) 

This representation of the Hermite polynomials is unfamiliar because of the operator 
produced by the composite one- variable similarity transformation between H^""^ and H^'^\ 
If this transformation is decomposed into elementary canonical transformations, direct 
evaluation leads to the more familiar Rodriques' formula for the Hermite polynomials. 
The decomposed transformation is 

^(»)(5) = -l-/ dqe'^'^e^'^ (ig'e-'«"«'g'". (87) 

ZTI J —oo J —oo 

This may be evaluated by first rewriting q'^ as (idg)'^ acting on the exponential exp{—iqq') . 
It may be extracted from the q' integral which then gives S{q). Integrating by parts n 
times transfers the idq operators to act on the remaining exponential terms 

^(«)(g) = r dq{-idqYe'^'i+^^'^5{q). (88) 

J —oo 

Completing the square of the argument of the exponential gives 

1 2 

— {q + 2iujqf+ujq , 

from which the purely q part can be extracted from the integral. The —idq derivatives 
act equivalently to —dq/2uj derivatives, and after converting them, they can be removed 
from the integral. This leaves a Gaussian integrated against a delta function which is 
immediately evaluated. The result is 

^/,(«)(g) =e"^'(-^)"e-'^'''. (89) 
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This is proportional to the Rodrigues' formula for the Hermite polynomials 

^^n(0 = e«'(-9g)"e-«^ (90) 

From this form, it is immediate that e^^ {—d^)e~^^ = —d^ + 2^ is the raising operator. 

There is a second linearly-independent solution of the harmonic oscillator which was 
not obtained by this canonical transformation. This solution is not normalizable, but 
from the standpoint of simply solving the differential equation, this is not important. A 
second canonical transformation which trivializes the super-Hamiltonian in a different 
way produces the other solution. That this should be necessary is not surprising since 
the original super-Hamiltonian is a second order differential operator while the trivialized 
one is only first order. 

Beginning from the harmonic oscillator Hamiltonian ([THD, the one- variable similarity 
transformation 

p = p"- - iuq"-, q = q"' 
is made to cancel the quadratic term in the coordinate, leaving 

The composite one- variable shift 

pa ^ q" = q'' - ip^/2uj 

cancels the quadratic term in the momentum, giving 

H^^'^ = -2iujqY - uj. (91) 
Finally, the point canonical transformation 

/ = e-«'p', = e"'. 

transforms this to action-angle form 

H^^^ = -2iup' - UJ. (92) 

The super-Hamiltonian can be trivialized by a two-variable similarity transformation, 
but this is unnecessary. 

The wavefunction ij)^^^ is given in terms of if)'^^'^ by 

^^°H9) = 5-c.,V2'^pV4-^ln,^^'^- (93) 

From the eigenf unctions of H^^'> 

^^1) = g-(«+l)'?-*(2n+lMo^ 
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one has 

These are clearly not normalizable for any n. 

The problem of the inverted harmonic oscillator with Hamiltonian 

can be solved by the above canonical transformations after uj is replaced by iuo. This is a 
scattering problem, so both independent solutions are delta-function normalizable, with 
no quantization of n. This emphasizes the importance of both canonical transformations 
to triviality. 



2.4 Time-Dependent Harmonic Oscillator 

The time-dependent harmonic oscillator can also be solved by canonical transformation. 
This has been done previously with a different sequence of canonical transformations by 
Brown[0. Here, the approach will be to parallel the solution of the time-independent 
harmonic oscillator in trivializing the super-Hamiltonian. This emphasizes the connection 
between the time- dependent and time-independent problems. It is conjectured that the 
parallel structure shown here carries over to time-dependent versions of other exactly 
soluble problems. 

One begins with the super-Hamiltonian 

n^^) =p^+p^ + uj\^ (96) 

where the angular frequency u = uj{qQ) is a function of time. The quadratic term in the 
coordinate is cancelled by making a two- variable similarity transformation generated by 
Fa = /(gQ)g^/2. This gives the super-Hamiltonian 

with the constraint 

if^^a + 2f = 2uj\ (98) 
This Ricatti equation is linearized by the substitution 

f = i-^ (99) 



which gives 



= -^uj'i,. (100) 
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To go further requires the specific time-dependence of u){qQ)- 

The next step is to cancel the quadratic term in the momentum with a second two- 
variable similarity transformation generated by = g{qQ)p°' /2. This gives the super- 
Hamiltonian 

T^W = + (2 - i^,^. + 4/^)//2 - 2ifqV - f (101) 

with the condition 

ig^^,-Afg^2. (102) 
This equation can be integrated to find 

g = exp(-i4 J fdq'o) J -2iexp(i4 J fdq',)dql (103) 

The super-Hamiltonian is then 

n^'^ = - 2ifqy - f. 

The coordinate change 

p& = e-«y, q'' = (104) 
eliminates the coordinate from the super-Hamiltonian 

n^<^)^pl-2ifp^-f. (105) 
Finally, the super-Hamiltonian is trivialized to 

(106) 

by the composite two-variable shift generated by 

F,= (-2p^ + i) j fdql. 
The wavefunction is given as 

^(1) ^ e-^^Pe^e-^^'e-^-V^^"^- (107) 
Since ip'^^'^ is any ^Q-independent function, this gives the result 

^(0) _ e/(.o)«V2eP(.o>V2p^^^g{-2p+0//(.o)'i«o^(i)(g). (108) 
The time-independent result is recovered when / = —uj and '0^^-*(g) = e"^. 
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2.5 Particle on 2n + 1-sphere 

As a final example, consider the radial Hamiltonian for a particle propagating on an 
2n + 1-dimensional sphere 

H^^^ ^p^ -2nicotqp. (109) 

Its solution by canonical transformations illustrates the use of the differential intertwining 
canonical transformation. First make the point canonical transformation 

= p, = cosq. 

sm q 

This is stated in inverse form, contrary to the convention followed above. This is often 
useful with point canonical transformations because it is more intuitive when looking for 
changes of variable to simplify an equation. The transformed Hamiltonian is found to be 

= (1 - q^\^^ + (2n + (110) 

This is recognized as the equation for the Gegenbauer polynomials. 
The differential intertwining transformation 

can now be used to cancel the n-dependence of the Hamiltonian. Noting that 

= i^^qVy - f + 2mq'-, - (HI) 
p p p° 

one finds 

i/W = (1 - g''')/' + igV - (112) 

Clearly, it would have been possible to shift n by any amount: this gives the relation 
between Gegenbauer polynomials of different n. 

Finally, undoing the original point canonical transformation 

p^^^^p, q^^cosq', 
sm q' 

gives the free-particle Hamiltonian 

(113) 

Because the physical problem was that of a free-particle on a sphere, this is the free- 
particle on a circle. The spectrum of H^^^ is discrete, and the constant shift produces 
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a time-dependent phase factor e*" '^o relative to the usual free-particle eigenfunctions on 
the circle 

= cosmge-^'"'''o. (114) 

(The other independent solution is found by using sinmg.) The original wavefunctions 
are given in terms of the free-particle eigenfunctions by 

V'^^ = Pcos,p"Parccos,^S„(g) = (^d\ cos(m + n)q e-(™^+2™„),o_ ^^^5) 

\smg ) 

The indexing is determined by the condition that m = corresponds to the normalizable 
solution with lowest energy. This agrees with the result obtained by the intertwining 
method^ and is recognized as a formula for the Gegenbauer polynomials c^^(cosg). 

In principle, this result is valid for real n. For n non-integer, one requires an inte- 
gral representation of the fractional differential operator. It is likely that this can be 
constructed by manipulating the definition of the composite similarity transformation 
in terms of the Fourier transform of an ordinary similarity transformation. This would 
be analogous to the discussion of the origin of the Rodrigues' formula for the Hermite 
polynomials. There are subtleties involving endpoints of the integrals which are beyond 
the scope of this paper. I hope to return to this in a later work. 



3 Conclusion 

It has been shown how, using a few elementary canonical transformations which have 
quantum mechanical implementations, a super-Hamiltonian can be trivialized and, thereby, 
its solutions found. The fact that the infinitesimal versions of these elementary transfor- 
mations classically generate the full canonical algebra is argued to imply that in principle 
any canonical transformation can be implemented quantum mechanically. Issues of op- 
erator ordering, to be sure, break the parallel structure between classical and quantum 
canonical transformations, so that in general different transformations are needed to 
reach the trivial super-Hamiltonian in each case. This raises the interesting possibility 
of the inequivalence of classical and quantum integrability. 

The method of intertwining was shown to be equivalent to constructing a canonical 
transformation between the two operators which are intertwined. This lays the foun- 
dation for explaining the relation between intertwining and Lie algebraic methods as a 
consequence of the fact that the dynamical symmetry group of a system is a subgroup 
of the infinite-dimensional group of canonical transformations. This explanation should 
extend as well to the infinite-dimensional symmetry groups behind the integrability of 
nonlinear systems like the Korteweg-deVries equation. These topics will be addressed in 
a later work ITT . 
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The key step in the construction of a differential intertwining operator was recognized 
as a particular composite similarity transformation, given the name of the differential in- 
tertwining transformation (^81) . This transformation has the remarkable property of 
being much simpler in the quantum case than in the classical, due to the noncommuta- 
tivity of the coordinate and momentum operators. Because of its role in the differential 
intertwining operator, this transformation is key to the integration of all differential equa- 
tions of hypergeometric type and their generalizations. More, the differential relations 
among hypergeometric functions of different indices and all formulae of the Rodrigues' 
type can be explained with it in terms of canonical transformations. 

The applications of canonical transformations discussed here only begin to explore 
the possibilities. My expectation is that the convenience of the language of canonical 
transformations will be found to be most powerful in problems in many-variables where 
conventional methods have been less well developed. This is a subject for future work. 
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P = -Q 
q ^ p' 



p ^ q 
q = -p' 



p = p'-if{q')^g, 
g = q' 



q = f{q') 



^ _ 1 ,f P,ip)¥'\q) = ^IZod<ie''''l-oodQe-''^^^^^^^ 



p = c/(p') 

Figure 1: Elementary and composite elementary canonical transformations 
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